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CN ■ We propose a systematic coarse-grained representation of block copolymers, whereby each block is 

reduced to a single "soft blob" and effective intra- as well as intermolecular interactions act between 
centres of mass of the blocks. The coarse-graining approach is applied to simple athermal lattice 
models of symmetric AB diblock copolymers, in particular to a Widom-Rowlinson-like model where 
1/^ ' blocks of the same species behave as ideal polymers {i.e. freely interpenetrate), while blocks of 

opposite species are mutually avoiding walks. This incompatibility drives microphase separation for 
copolymer solutions in the semi-dilute regime. An appropriate, consistent inversion procedure is 
used to extract effective inter- and intramolecular potentials from Monte Carlo results for the pair 
' distribution functions of the block centres of mass in the infinite dilution limit. 
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^ I. INTRODUCTION 

^ ' Recent years have witnessed considerable theoretical interest in statistical mechanics of multicomponent systems 
O I involving multiple length and time scales. To be tractable, statistical descriptions of such multi-scale systems must 
resort to controlled coarse- graining methods. One widely used coarse- graining strategy for complex fluids is to 
determine effective interactions between large ("dressed") particles {e.g. macromolecules or colloidal particles) by 
systematically tracing out the degrees of freedom of small particles {e.g. counter- ions or solvent molecules 
^ ' A particularly successful application has been to dilute and semi-dilute polymer solutions: effective interactions 
, between polymer centres-of-mass (CM) are determined by takiiig appropriate averages over monomer conformations 
' for fixed distances r between the CMs of interacting linear0|[3|a]|3 or star polymersQQ. In the case of linear 
J — ' polymers, to which the subsequent discussion will be restricted, the resulting CM — CM pair potential for polymers 
in good solvent, modelled by self avoiding walks (SAW) on a cubic lattice, is very soft, repulsive and finite at all 
lO ] separations r, with a range of the order of the radius of gyration Rg of the polymer coils. The general shape is 
I reasonably well represented by a Gaussian of amplitude on the order of fc^T, reflecting the essentially entropic nature 
of the effective interaction. While the overall form of the effective pair potential is not very sensitive to polymer 
concentration under good solvent or high temperature conditions, it turns out to be very sensitive to concentration 
and temperature when the latter is lowered toward the 9 solvent regime 6J 9]. 

Replacing a full monomer level representation of polymer coils by a "soft particle" description based solely on 
effective pair potentials between CMs represents a considerable reduction in the number of degrees of freedom, and 
Q ' hence leads to a concomitant reduction in computational effort in simulations of large-scale phenomena involving 
O ' many polymers. Moreover, use of spherically symmetric effective pair potentials allows a direct exploitation of the 
theoretical arsenal developed over the years for the study of the bulk|^3| and interfacial^J properties of simple 
• ^ fluids. Last, but not least, inspection of the effective pair potentials leads to new insights into the phase behaviour^ 
' and stabilityl^ of macromolecular solutions. 

In this paper we extend the above coarse-graining strategy to the case of solutions of block copolymers and more 
speciflcally to symmetric diblock copolymers made up of two strands A and B of equal length. AB copolymers will 
be reduced to "soft diatomics", with inter- and intra- molecular interactions between sites associated with the CMs of 
A and B strands. The nature of A-A, A-B and B-B effective pair interactions depends on solvent conditions. Most of 
the results presented in this paper will be based on a simplified model inspired by the celebrated Widom-Rowlinson 
(WR) model for fluid/fluid phase separation T^. In the macromolecular extension of the WR model, A and B strands 
are mutually avoiding {i.e. monomers of opposite species cannot overlap) while A and B strands will separately 
behave like ideal (Gaussian) chains, i.e. they will freely interpenetrate same species strands on different copolymers. 
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II. MODEL AND COARSE-GRAINING STRATEGY 



We consider lattice models of polymers, where monomers occupy the sites of a periodically repeated simple cubic 
lattice of size A^; the bond vectors connecting successive monomers point along the x, y or z directions only. Each 
polymer comprises A4 monomers and hence L = A4 — I bonds (or segments), such that the length of the polymer 
is L 6, if 6 is the segment length. A diblock copolymer AB is made up of two strands comprising Ma and Mb 
monomers respectively, with the M^th monomer of strand A connected to the first monomer of strand B. The total 
length is {La + Lb + l)b = {Ma + Mb — l)b. Restriction will be made here to symmetric copolymers, such that 
Ma — Mb = M and La = Lb = L. We will consider N such copolymers on an A'^ lattice, such that the A and B 
monomer concentrations are ca — cb — MN/A^ and the overall monomer concentration is c = ca + cb- Individual 
copolymer chains are characterised by three radii of gyration RgA, RgB and Rg, defined as usual by 



Ma 

Ma 



i—1 

where l^f is the position of the i*'' monomer of species A on the lattice, while Ra = X^i^i '^t I^A is the position of 
the centre-of-mass {CM) of strand A. The statistical average < . > is taken over properly weighed chain conformations. 
Similar definitions hold for RgB and for the overall radius of gyration Rg] the overall CM is R = {Ra + Rb)/'^- 
Another important length is the root mean square distance between the CMs of the two strands of a copolymer: 



An elementary calculation shows that: 



Rab = [<\Ra - Re? >Y'^ (2) 



^ 2^1a+ ^R^gB + ^R\b- (3) 



The A-A, A-B and B-B interactions between non-adjacent monomers of the same or opposite species determine the 
structure and phase behaviour of individual copolymers and of many-copolymer systems. A number of models have 
been examined in the polymer literature: the more interesting are those which account for some degree of incompatibil- 
ity between A and B monomers, which leads to microphase separation 13J 14] . In this paper we focus on the simplest, 
athermal models combining ideal (or random walk, with backward overlapping steps allowed) strands (I) and excluded 
volume (self or mutually avoiding walks) strands (S). The eight possible combinations of monomer interactions will 
be labelled by three indices JKL, with each index either I or S. The first index refers to the interaction between A 
monomers (on the same or different polymers), the second index refers to the A-B cross interaction, while the last 
index refers to the B-B interaction. Of the eight possible combinations, two (III and SSS) reduce to homopolymers 
of molecular weight 2Af, while IIS and SIT as well as ISS and SSI are degenerate for the symmetric {Ma = Mb) 
copolymers considered here, leaving four distinct, non-trivial copolymer models, namely IIS, ISI, ISS and SIS. Among 
these, the ISI model may be regarded as the macromolecular equivalent of the WR model J^J. The macroscopic phase 
separation occurring in the latter is prevented here because A and B strands are tethered, but microphase separation 
into a lamellar structure may be expected. Most of the subsequent discussion is hence devoted to the ISI model. 
Physically, self avoiding walks and random walks (ideal polymers) reflect solvent conditions. The former model corre- 
sponds to good solvent conditions associated with swollen chains, while random walk polymers are generally considered 
as approximating 0-solvent conditions, at least for low polymer concentrations 15j. Thus in the ISI model, A and B 
strands behave individually as if they were in a 0-solvent, where monomer excluded volume effects are compensated by 
nearest-neighbour, solvent-induced attraction, while A-B interactions are purely of excluded volume nature (mutually 
avoiding walks). A more general model than those considered here contains both excluded volume interactions (single 
occupancy constraint) and nearest neighbour interactions ea(3{a,(3 — A or B), thus introducing an explicit tempera- 
ture dependence. This standard model has been extensively studied both in the low concentration(single copolymer) 
limit(16i|. and in the melt JJ:]. The present investigation covers special cases of the generic model for finite copolymer 
concentrations, in the dilute and semi-dilute regimes. The key characteristic of these models, and in particular of the 
ISI model, is their extreme non-additivity, as in the WR model. 

The range of concentrations which will be considered is best characterised by the ratio of the copolymer density, 
p — N/{Ab)^ over the overlap density p* = 3/(47ri?g'^), where Rg is the radius of gyration at zero density. Dilute 
solutions correspond to p/p* <c; 1, while the semi-dilute regime is for p/p* > 1. 

The coarse-graining strategy which we shall pursue generalises that developed earlier for linear homopolymers0|0@. 
An AB block copolymer will be represented by two spherical "blobs" centred on the CMs Ha and Rb, and tethered 
by an anharmonic entropic spring deriving from an intramolecular potential (I)ab{'''), where r = \Ra ~ R-sI- The 
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A and B blobs on different copolymers will interact via effective pair potentials VAA{f),VAB{f) and VBsif) acting 
between the CMs. As in the earlier work on homopolymers, these 4 effective potentials are determined by an inversion 
procedure based on the CM — CM pair distribution functions SABir) (intramolecular) and gAAir), gAB{r) and gBB{r) 
(intermolecular). The latter are determined by Monte Carlo (MC) simulations of the initial, monomer- level model. 
In this paper these pair distribution functions and resulting effective pair potentials will be determined only in the 
low concentration limit, by simulating copolymer pairs. While the inversion procedure is trivial in that limit in the 
case of homopolymers (represented by single blobs), it turns out to be much more involved in the copolymer case, 
since four effective blobs are involved, as will be discussed in more detail in section IV. 

III. MONTE CARLO RESULTS 

With the above coarse-graining objective in mind, we have carried out extensive MC simulations of symmetric AB 
diblock copolymers on a cubic lattice. Most simulations were for the ISI model, but some preliminary results have 
also been obtained for the IIS, ISS and SIS models; test runs were carried out for the SSS model to compare with 
previous results of the equivalent SAW homopolymers. 

MC runs were carried out on periodic lattices of size at least A = 100. The numbers of copolymers were adjusted 
to achieve any given polymer density p/ p*] thus N was taken to be 1 to determine infinite dilution, single polymer 
properties, iV = 2 if effective pair interactions were to be determined in that limit, while for p/ p* > 1, N varied 
between a few hundred and over a thousand. The number of monomers 2M in each strand of a copolymer was varied 
between 60 and 1000. Configuration space was sampled by using pivot and translation moves [T^|l8||. polymers were 
subjected to between 5 x 10^ and 10^ MC moves, depending on the total number of polymers in the system. 
Results for the 3 radii of gyration and the mean CM — CM distance Rab are collected in Table 1, for the four 



2M=320 



Type 


Rg 


RgA 


RgB 


Rab 




ISI 


8.19 


5.27 


5.28 


14.11 




IIS 


9.99 


5.16 


8.53 


12.47 




ISS 


10.85 


5.27 


8.69 


16.08 




SIS 


12.13 


8.55 


8.50 


17.11 




SSS 


12.96 


8.70 


8.64 


19.14 




2M=500 


Type 


Rg 


RgA 


RgB 


Rab 


B2 


ISI 


10.228 


6.588 


6.586 


18.23 


4.84 


IIS 


12.902 


6.455 


11.146 


15.60 


3.76 


ISS 


13.872 


6.560 


11.279 


20.66 


5.71 


SIS 


15.798 


11.148 


11.146 


22.34 


3.87 


SSS 


16.831 


11.273 


11.271 


24.96 


5.70 


2M=2000 


Type 


Rg 


RgA 


RgB 


Rab 




ISI 


20.47 


13.20 


13.18 


40.32 




IIS 


28.94 


13.13 


25.54 


31.28 




ISS 


30.33 


13.13 


25.59 


45.23 




SIS 


35.95 


25.38 


25.36 


50.44 




SSS 


38.39 


26.06 


26.04 


56.23 





TABLE I: Radii of gyration Rg, mean intramolecular distance Rab and second virial coefficient B2 in units of Rg (2M=500 
only) for ideal/SAW systems. 

copolymer models as well as for the SSS model, which serves as a check against the well documented SAW homopolymer 
model, for which Rg ~ 0.44L'' with v ~ 0.587 in the scaling limit To*!. The trends in the various quantities in Table 1 
for the various models agree qualitatively with expectation. The internal CM-CM pair distribution function sab and 
the resulting "entropic spring" potential i. e. the effective bonding potential between the CMs of the A and B strands 
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of a single copolymer: 

(l)ABir) = -ksTlnsABir) (4) 
are plotted in figure While the distribution functions SAsir) for the IIS and SIS models peak at the origin, as one 




gm gm 



FIG. 1: left: Zero density limit of the intramolecular CAf-pair distribution function SAsir) for different models, versus r/Rgm, 
where Rgm ~ {RgA + RgB)/'2 is the mean radius of gyration of the A and B strands, right: The corresponding intramolecular 
pair potentials (pAsir), in units of fcsT, defined in equation|l] 

might expect since both strands can freely interpenetrate in both cases, the SAsir) have peaks at r > Rgm, where 
RgTu = {RgA + RgB)/'^, in the ISI, ISS and SSS models, because of the single occupancy constraint on monomers of 
opposite species in all three cases. These three copolymer models behave hence like macromolecular diatomics, i.e. 
they are, on average, elongated, as confirmed by inspection of the eigenvalues of their tensor of inertia (not shown 
here). The corresponding intramolecular potentials c/jab exhibit a maximum at full overlap (r = 0) and a minimum 
for r ~ l.bRgrn. Beyond that minimum, all potentials increase rapidly, illustrating the "entropic spring" character 
of the effective interaction between the CMs of the two strands. Because the two strands are tethered, (t)AB must 
diverge when r — > L x 6, corresponding to the highly improbably fully stretched conformation of the AB copolymers. 
Note that relative to the r = value, the potential minimum is considerably deeper for the ISI model, compared with 
the ISS and SSS models. In other words elongated conformations have higher probability for the ISI model. This is 
because the two ideal strands in the latter model are more compact than self avoiding strands, as reflected in their 
shorter radii of gyration. The potential barrier to interpenetration of mutually avoiding strands is thus larger when 
the strands are individually ideal, since the probability of the forbidden overlap of monomers of opposite species is 
larger. The results discussed so far are for isolated copolymers. We henceforth restrict all further considerations of 
many polymer systems to the ISI model. 

We first investigate the infiuence of polymer concentration on the intramolecular CM distribution function SAB{f). 
Figure 121 shows the MC data for the ISI model at four different concentrations, as well as the corresponding Rab, 
which is related to the properly normalised sab{t) by: 

R\b = 47r / SAB{r)r^dr. (5) 

SAB{r) is seen to be rather insensitive to concentration in the dilute regime as one might intuitively expect, Rab 
contracts slightly as p/p* increases before expanding at higher densities where microphase separation sets in. In terms 
of the coarse-grained diatomic analogy, this means that the two-blob "molecule" hardly contracts upon compression, 
and that this coarse-grained entity remains a valid concept in the semi-dilute regime. Note however that the "entropic 
spring" potential 4>ab is no longer given by eq^at finite concentration, because the problem ceases to be a purely 



5 




FIG. 2: Intramolecular CM pair distribution function SAs(r) of the ISI model as a function of r /RgA at four densities p/ p*. 
The inset shows Rab (defined in equation|SJ versus p/ p*. 



two-body one. A more complex inversion procedure is required, to which we return in section IV. 
In order to gain more insight into the physical significance of the ISI model, we have determined the osmotic equation 
of state of the model in the dilute and semi-dilute regimes. The equation of state can be efficiently and accurately 
determined in a single MC run, by subjecting the copolymers to a gravitational field until sedimentation equilibrium is 
reached, determining the resulting inhomogcncous monomer or CM concentration profile p{z) (where z is the altitude) 
and extracting the osmotic pressure P from the measured profile by invoking hydrostatic equilibrium pl f po | 

dP(z) 

— — = -mgp{z) (6) 
az 

where m is a fictitious buoyant mass and g the acceleration due to gravity; these parameters may be adjusted in the 
simulation, such that the resulting sedimentation length C = ksT /mg is larger than the polymer radius of gyration 
Rg , thus ensuring the validity of the macroscopic equation |^ on the mesoscopic scale set by i?g . A snapshot of a 
typical configuration in the sedimentation column is shown in figure 13 while the concentration profile p{z) averaged 
over millions of such configurations generated in the MC run is plotted adjacent in the figure. The configuration 
in figure |21 clearly hints at the existence of microphase separation of the ISI model for copolymer concentrations 
pI P* ^ 2. Horizontal cuts through the sedimentation column, shown in figureQ] exhibit alternatin g stripe s of A and 
B strands, characteristic of the lamellar phase observed in many symmetric block copolymer melts p^[l4|. However, 
while the control parameter in melts is the temperature, the control parameter for the athermal ISI model is the 
polymer concentration. 

Further evidence for the microphase separation comes from the osmotic equation of state, extracted from the 
concentration profile as explained above. The dimensionless equation of state Z = P/ pksT, computed for copolymers 
with total number of monomers 2M=120, 250, 500, and 2000 are plotted versus p/ p* in figure [SJ no significant size 
dependence of Z is seen. Z increases from its infinite dilution limit 1 up to p/ p* ~ 2.5, where it reaches a maximum, 
before decreasing slowly to what appears to be an asymptotic value for p/ p* > 10. With increasing density the system 
moves from a disordered region to one where the lamellar order increases, and there are fewer interactions between A 
and B strands. The excess pressure(compared to ideal polymers) is decreases in this region, as a consequence of the 
lamellar ordering minimises the number of A-B overlaps. The initial slope of Z{p) agrees well with the second virial 
coefficient reported in table 1, while at high concentrations, the osmotic pressure increases sub-linearly with p/ p*. 
In future work the anisotropic pair structure will be calculated to characterise the onset of microphase separation in 
more detail. 
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FIG. 3: Left: typical configuration of N=1600 ISI copolymers under gravity {(^ = kBT/mg — 1.66-Rg). The copolymers are 
confined to a vertical box, open-ended in the z- direct ion (2 > 0), and periodically repeated in the horizontal x and y directions, 
the horizontal base is 100x100 lattice units. A and B strands are pictured as black and grey, spheres of radius R = RgA- 
Microphase separation into a striped pattern in clearly visible at lower altitudes. Right: monomer density profile p{z) for the 
system, as shown left. 



IV. SOLVING THE INVERSE PROBLEM 



We now return to the main objective of the present paper, namely to determine the efFective pair potentials VAA{r), 
VAB{r) and VBsir) between the CMs of A and B strands on different copolymers. In the case of the symmetric ISI 
model, VAAif) = vbb{^)- From earlier work on homopoIymersQ, which have coarse-grained interaction characterised 
by a single effective pair potential, we expect the effective pair potentials to be state dependent i.e. to depend on 
polymer concentration. The subsequent discussion will be restricted to the limit of vanishing concentration, i.e. to an 
isolated pair of copolymers. To this end the intramolecular pair distribution functions gAAir) = gBsir) and gAsif) 
have been computed by MC simulations of a pair of AB copolymers. The procedure is similar to that used previously 
for homopolymers'5] G], whereby conformations of the two isolated polymer coils are generated, and the CMs of the 
two coils are moved toward each other, while histograms are collected of the allowed configurations as a function of the 
CAI — CM distance. In the case of diblock copolymers, care must be taken to correctly generate the partial gap{r). 
The simulation must be adjusted to move polymers such that the relevant CM^ — CMp distance is changed. The MC 
results for the three low density pair distribution functions of the ISI model are shown in figure together with the 
pair distribution functions of the overall CMs of the copolymers, gcci'i')- These results together with the previously 
determined intramolecular pair distribution function sab(''); are then inverted to yield effective, monomer-averaged 
pair potentials VAAir) — VBsir), VABir) and vccif)- While the inversion oi gccii^) to yield vcci'i') is trivial; 



vcc{r) = -ksTlngccir) 



(7) 
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FIG. 4: Two horizontal cuts of width = RgA of the configuration in figure 13 at altitudes z = 5.5-Rg(left) and z — 
9.5-Rg(right). A and B strands are again pictured as black and grey spheres of radius Rg. The lamellar structure is starting to 
break up at 2 = 9.5Rg. 



o^2M=120 




2 4 6 8 10 

p/p* 

FIG. 5: Osmotic equation of state Z = P/{pkBT) of the ISI model as a function of p/p* as extracted from equation|^ The 
results are for copolymers of different lengths 2M=120, 250, 500 and 2000 



the inversion of gAAir) and gAsir) to yield VAAir), VABii"), given SAsir), is considerably more complicated, since 
one is in effect facing a four-body problem. For site-site interaction models of "diatomics" , the exact relation between 
gapij) and Vapir), for a given SAsif) is in the low density limit : 
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\im hAA{r) = fAA{r) + [1 + .fAA{r)] '^^ [/ab(x)sba(x - r) + sab(x)/ba(x - r)] 

+ y"c?x J dysAB{^)sBA{y - r) [/bb(x - y) + fAB{y)fBB{^ - y) + /sa(x - r)/BB(x - y) 

+/As(y)/sA(x - r) + /As(y)/i3A(x - r)/BB(x - y)] |> (8) 



\hnhAB{r) = fAB{r) + [1 + | y [/aa(x)sab(x - r) + sab(x)/bb(x - r)] 

+ y c?x y dysAB(x)s^s(y - r) [/ba(x - y) + /AA(y)/BA(x - y) + /bb(x - r)/B^(x - y) 

+/AA(y)/i3s(x - r) + fAA{y).fBB{^ - r)/sA(x - y)] ^ . (9) 



where /i^^ (r) = (r) — 1 and the /^^ (r) are the Mayer functions 

/a/3(r) = /;3a(r) = exp [-Pvapir]] - 1, (10) 

The inversion procedure now amounts to solving the coupled integral equations El and |51 for the fafsir), using the 
haf){r) and Sapir) from the MC simulations as input. Once the fapir) have been calculated by an appropriate iterative 
solution of the two coupled integral equations, the effective intermolecular potentials Vaf}{r) follow from eauation llOl 
The numerical solution of equations |H1 and is facilitated by the fact that all integrals on the r.h.s. are convolution 
integrals, and hence easily evaluated by Fourier transformation, except the last term of the double integrals, which 
involves five factors and corresponds to a fully connected "bridge" diagram, which cannot be resolved by Fourier 
transformation. If the bridge term is left out, in the spirit of the familiar hyper-netted chain (HNC) approximation[Toj . 
major difficulties are encountered when attempting a numerical solution, because of the inconsistency introduced by 
neglecting the bridge term, which is a violation of the connectivity constraints. However the task of including the 
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FIG. 7: Optimal representation of the MC data for A-nr^hAA^r) and 47rr^/iAs(r) of the ISI model in the zero density limit, 
achieved by the inversion procedure of equations I8I9I and IIUI The unknown fapir) are parametrised by sums of 10 Gaussian 
functions. 




FIG. 8: Left: Mayer functions fAAir) — fBB{r) and fAB{r) versus r/Rgm, as derived by the inversion procedure outlined in 
the text and in figure|7| Right: Resulting effective pair potentials VAA{r) = VBB{r), VAB{r) and vccir), versus r/Rgm', the 
effective potentials correspond to the zero density limit. 



bridge term may be achieved by noting that all integrals in equations|Sland|51may be calculated analytically, if SAsi^), 
fAAir), fAsir) are sums of Gaussian functions centred on r = 0. We have hence fitted the MC data for SABir), shown 
in figure 121 by a sum of four Gaussian functions. The (two) unknown Mayer functions fAAir) = fsBir), fAsir) are 
represented by a sum of 10 Gaussian functions, and the amplitudes and widths are varied until the resulting pair 
correlation functions Anr'^hapir) yield the best fits to the MC data as illustrated in Fig[7| The corresponding optimal 
Mayer functions fAAir), fAsir) and the resulting effective pair potentials VAAir) and VAsir) are plotted in Fig|Sl 



10 



As expected, VAA{r) = Vbb{t) is small compared to fc^T for all site-site distances, while VABir) is roughly Gaussian 
in shape, with an overlap (i'abC'" = 0)) value of i.iksT. The effective pair potential vccif) between the overall 
copolymer CMs, calculated from eq. Eland the MC data for gcc{f), is also shown in FiglHl Its amplitude at full 
overlap, vcci''' = 0) ~ 1.25fcBT is almost three times smaller than VABif — 0) but its range is significantly wider, as 
one may expect, since the effective potential vcci''') involves an implicit averaging over relative orientations of the 
elongated diblock copolymers. It is interesting to note that the entropic barrier of 3.2kBT for complete overlap of the 
mutually avoiding strands A and B of the two copolymers is substantially higher than found for homopolymers, where 
v{r = 0) « l.SfcsT for sufficiently long polymers^J^Sj 9]. This is a consequence of the fact that monomers of the same 
strand can overlap (ideal polymers), and that the strands are hence more compact than the swollen self-avoiding walk 
polymers. Interestingly, the effective potential ^^^(r) for two copolymers is rather close to that found for binary 
mixtures of untethered A and B polymers 12^. 



V. CONCLUSION 



We have introduced and investigated a highly simplified model of a symmetric diblock copolymer which leads to 
microphase separation, the ISI model. The binary mixture counterpart of ISI, where A and B strands are untethered 
is the macromolecular equivalent of the Widom-Rowlinson modeljlj, and like the latter leads to macroscopic phase 
separation 23). Despite the simplicity of the ISI model, a fully quantitative investigation of its structure, thermody- 
namics and phase behaviour would still be a daunting task, and only partial results have been reported in the present 
paper. In order to go further we propose a coarse-graining strategy, whereby an AB block copolymer is schematized by 
two "blobs" tethered by an entropic spring, deriving from an effective intramolecular pair potential ipABir) between 
the CMs of the two strands. We were able to extract (t)AB{r) and the intermolecular potentials VAAir) and VABir) 
from the site-site pair distribution functions Syis(r), gAAir), gABir) in the zero density limit, as generated by MC 
simulations for an isolated pair of interacting copolymers with 2M monomers (10^ ^ M ^ 10'^). The simulation of 
the resulting coarse-grained "soft diatomic" representation of the ISI model is orders of magnitude faster, since the 
number of degrees of freedom is reduced by a factor M. Preliminary investigations show that the coarse-grained model 
indeed leads to phase behaviour very similar to that of the original ISI model. Note that, by construction, the "soft 
diatomic" model will lead back to the exact pair distribution functions of the ISI, at least at zero density. From our 
earlier experience with effective pair potentials between the CMs of homopolymers, we expect some concentration 
dependence of the effective potentials for AB copolymers. In the case of the former, the inversion of the single pair 
distribution function at finite polymer concentration is straightforward, via the HNC closure which is extremely ac- 
curate for soft potentials'],^]. A similar inversion scheme for diblock copolymers could in principle be based on the 
"reference interaction site model" (RISM) formalism[^, with HNC closure, but the inconsistencies of this formalism 
are well documented, and will require careful attention for its implementation in an inversion procedure. Since no 
dramatic density-dependence of the effective potentials is expected, at least for the athermal ISI model, an excellent 
first approximation will be to use the zero density potentials derived in the present paper to investigate the behaviour 
of dilute and semi-dilute copolymer solutions. A severe check will be provided by a direct comparison of the pair 
distribution functions obtained at finite densities with the effective "soft diatomic" model, and from the full underly- 
ing ISI model. Work along these lines is in progress. The extension of the present work to asymmetric {Ma Mb) 
diblock, triblock and more complex copolymer systems is, at least in principle, straightforward. 
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